On the 2D behavior of 3D MHD with a strong guiding field 
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The Magneto-hydrodynamic (MHD) equations in the presence of a guiding magnetic field are 
investigated by means of direct numerical simulations. The basis of the investigation consists of 9 
runs forced at the small scales. The results demonstrate that for a large enough uniform magnetic 
field the large scale flow behaves as a two dimensional (non-MHD) fluid exhibiting an inverse cas- 
cade of energy in the direction perpendicular to the magnetic field, while the small scales behave 
like a three dimensional MHD-fluid cascading the energy forwards. The amplitude of the inverse 
cascade is sensitive to the magnetic field amplitude, the domain size, the forcing mechanism, and the 
forcing scale. All these dependencies are demonstrated by the varying parameters of simulations. 
Furthermore, in the case that the system is forced anisotropically in the small parallel scales an 
inverse cascade in the parallel direction is observed that is feeding the 2D modes k« = 0. 
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I. INTRODUCTION 

The existence of magnetic fields is known in many 
astrophysical objects, such as the interstellar medium, 
galaxies, accretion disks, star and planet interiors and the 
solar wind [J. In most of these systems, the magnetic 
fields are strong enough to play a significant dynami- 
cal role. The kinetic and magnetic Reynolds numbers 
involved in these astrophysical bodies are large enough 
so that the flows exhibit a turbulent behavior with a 
large continuous range of excited scales, from the largest 
where energy is injected, toward the finest where energy 
is dissipated. In many cases, strong large-scale magnetic 
fields are present that induce dynamic anisotropy in the 
small scales. Direct numerical simulations that examine 
in detail both large and small scale turbulent processes 
in astrophysical plasmas are very difficult to achieve, and 
only modest scale separation can be reached even with 
today's super-computers. One way around this difficulty 
is to model the effect of the large scale field by a uni- 
form magnetic field Bq, and thus study the small scales 
separately. 

In the presence of a strong uniform magnetic field 
the evolution of the turbulent fluctuating fields can be 
treated within the framework of weak turbulence theory 
(WTT). In this approach the nonlinearities are treated 
perturbatively, resulting in a slowly varying amplitude of 
the linear wave modes that in this case are the Alfven 
waves supported by the uniform magnetic field [U, 
However, the limit Bq — > oo is non-trivial because dif- 
ferent limiting procedures can lead to different results. 
Thus in order for the results of weak turbulence theory 
to hold a domain of size L (in direction of the field) 
sufficiently large needs to be considered, so that modes 
with small wave numbers in the direction of the field that 
satisfy the "quasi-rcsonancc" conditions are present. In 
particular WTT is valid when the following condition is 
met: 



(see ||) where \uk\ is the amplitude of the velocity field 
with wavenumbers ~ k (with fc|| and k± the parallel and 
perpendicular projection of k to the magnetic field re- 
spectively). L|| is the domain size in the direction of the 
field. The inequality on the right implies sufficiently weak 
nonlincarity while the inequality on the left is needed for 
the presence of quasi resonances. In this case WTT pre- 
dicts that the energy spectrum is proportional to kj 2 . 
On the other hand if 
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the system becomes "slaved" to the 2D modes k« 
that evolve independently [|| . Thus if the limit Bq — > oo 
is taken keeping the domain size L fixed the system be- 
comes two dimensional 0, Note that the two con- 
ditions ([U,© allow the existence of an intermediate 
regime. 

Here the large Bq limit is explored further by 
means of numerical simulations. Numerically, magneto- 
hydrodynamic (MHD) turbulence has been investigated 
by various groups in the last decade |6-10|. The results 
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of WTT were first demonstrated m 
transition to two dimensional dynamics has been inves- 
tigated more recently in [l3[. In all these investigations 
the flow was forced in the largest scale of the system. 
In this work the case where the system is forced in the 
small scales is explored and the possible development of 
an inverse cascade is examined. 

An inverse energy cascade is known to exist in two- 
dimensional (2D) hydrodynamic turbulence (T34l6| . as a 
consequence of the conservation of vorticity. It results 
in a /c~ 5 / 3 spectrum for the large scales and a k~ 3 sp ec- 
trum for the small scales. Strongly rotating flows p7H2l| 
or flows in thin boxes [22[ have been shown to have a 
dual cascade of energy with the large scale flow behaving 
like 2D with an inverse cascade while the small scales 
being three dimensional with a direct cascade. How- 
ever, 2D MHD turbulence does not conserve the vorticity 
and energy is cascading to small scales. On the contrary 
the square of the vector potential is cascading to larger 
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scales [23|, |24|. For the system under investigation wc 
cannot a priori predict if the flow under the influence of 
a strong magnetic field will act as a 2D- hydrodynamic 
flow by suppressing all magnetic fluctuations and thus 
have an inverse cascade; or if magnetic fluctuations per- 
sist and the system will act as a 2D-MHD flow and thus 
not exhibit an inverse cascade of energy. It is noted that 
a strong magnetic field is used very often to make flows 
of liquid metals behaving like two dimensional flows in 
experiments where an inverse cascade has been observed 
[2a . These flows however have very small magnetic 
Reynolds numbers and magnetic fluctuations are strongly 
suppressed. 



RUNS 


1 

G K /W 


v A 


k f L 


M 


Isotropy 


Rl 


5.0 


5.0 


8-10 


0.0 


I 


R2 


2.5 


2.0 


8-10 


0.0 




R3 


2.5 


10.0 


8-10 


0.0 




R4 


10.0 


5.0 


4-5 


0.0 




R5 


2.5 


5.0 


16-20 


0.0 




R6 


3.3 


5.0 


8-10 


0.4 




R7 


2.5 


5.0 


8-10 


0.6 




R8 


2.5 


5.0 


8-10 


0.0 


F(k||,0) = 


R9 


16.0 


5.0 


8-10 


0.0 


F(0,k ± ) = 



TABLE I. Table with the parameters of all runs. "I" in the 
last column stands for isotropic forcing. With boldface are 
marked the parameters that are varied with respect to Rl. 



II. FORMULATION AND NUMERICAL 
SIMULATIONS 

We consider a flow of a conducting fluid inside a triple- 
periodic box of size 27rL in the presence of a strong guid- 
ing magnetic field Bq in the z direction. The system is 
forced by a mechanical force f and an electro- motive force 
£ . The non-dimensional MHD equations then read: 

a t u + u-Vu= V A d z b + b • Vb - VP + G~J V 2 u + F 

d t h + u • Vb = V A d z u + b • Vu + G~J V 2 b + MV x E, 

where u is the velocity field and b is the magnetic field. 
Both fields are assumed to be solcnoidal V • u = V • b = 
0. F = f/||f|| is the external mechanical force normal- 
ized to unit amplitude (where || • || stands for the Li 
norm) . E is the external electro-motive force normalized 
so that its curl has unit amplitude E = £/||V x £\\L. 
The equations have been non-dimensionalized using the 
box size L and the forcing amplitude ||f||. With this 
choice four non-dimensional parameters appear. G K is 
the kinetic Grasshof number G K = ||f||L 3 /j/ 2 , where v 
is the viscosity. G M is the magnetic Grasshof number 
G M = ||f||L 3 /77 2 , where 77 is the magnetic diffusivity. In 
all the runs performed in this work G K = G M . The 
amplitude of the external magnetic field relative to the 
forcing is controlled by the parameter V A = Bq/ y/\\f\\L. 
Finally M = [V x £||/||f|| expresses the ratio of the 
electro-motive to the mechanical forcing. 

A possible alternative to this non-dimensionalization 
choice would be to use the kinetic and magnetic Reynolds 
numbers, typically defined as Re = ||u||L/z/ and Rm = 
\\u\\L/r) respectively For our problem at hand however 
where an inverse cascade is present it is not an attractive 
choice because the amplitude of the velocity is changing 
with time throughout the duration of the computation 
rendering Re a function of time. 

Both forcing mechanisms used in the numerical simu- 
lations consisted of a sum of Fourier modes with wave- 
numbers inside a spherical shell |k| = kf. The phases 
of the modes were changed randomly every time interval 
r ~ y/L/\\f |. The forcing is isotropic in all runs except 
the last two that we discus in section llil El There was no 



averaged helicity or cross-helicity injected by the forcing 
by choosing: (F • V x F) = and (F • V x E) = 0. Since 
we are interested in the presence of an inverse cascade we 
are forcing relatively large wave numbers. 

The MHD equations were solved using a standard 
pseudo-spectral method and a third order in time Runge- 
Kuta [13, [28J. The resolution used in all the runs was 
512 3 grid points. In the 9 different runs that were per- 
formed the amplitude of the external magnetic field, and 
the way the system is forced was varied. The table HI 
gives all the parameters of the runs. 

The choice of G K in most runs is rather conservative 
because it is not known beforehand what effect on the 
resolution requirements a change in each of the parame- 
ters has. In any case, in all runs a well resolved spectrum 
was observed. The last run (R9) has a large value of G K 
because in this case for similar forcing amplitude with 
the other runs the flow is less efficient in absorbing en- 
ergy because the k\\ = modes are not forced. In this run 
also the time scale of the forcing was set to r ~ (£?ofc/) _1 
to be closer in resonance with the forced Alfven-modcs 
and to improve this absorption efficiency. 

The diagnostics used are based on energy spectra and 
energy fluxes that are now defined. If Uk and bk are 
the Fourier modes of the velocity and magnetic field of 
wavenumber k then the two dimensional energy spectra 
E u (kj_, fen) and Eb(kj_, fen) are defined as 

where the sum is restricted in the wavenumbers fc|| < 
\k z \ < fc||+l and k± < J~k%. + ky < k± + l. The averaged 

energy spectra in the parallel direction are then defined 
as 

E u (k x )=J2 E ^' k \0> E b (k ± ) = J2E b {k ± ,k l{ ). 

fc,i fen 

The total kinetic energy E K and magnetic energy E M are 
then given by 

E K = E u (k ± ), E M = E b (k±). 
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FIG. 1. Time evolution of the kinetic energy E K (solid line) FIG. 2. The energy flux for run Rl in the perpendicular 
and the magnetic energy E M (dashed line) for run Rl. direction (solid line) and parallel direction (dashed line). 



Since the development of a direct or an inverse cascade 
is expected to be anisotropic we need to define the flux of 
energy through an arbitrary surface in Fourier space. If 
u D and b D stand for the projection of the two fields u, b 
to the flows whose Fourier transform contain modes only 
inside the Fourier domain "D" , the energy flux through 
this domain is given by 

U D = J[u D u ■ Vu-u D b • Vb+b D u • Vb-b D b • Vu]dV. 

See [29| for more details. In this work we will consider 
the flux through cylinders and planes. By II j_ (k± ) we will 
refer to the flux of energy through a cylinder of radius 

\Jk^, J r ky = k± and by n_i_(fc|| ) we will refer to the flux of 

energy through the planes \k z \ = fen. Positive flux implies 
cascade of energy to the small scales while a negative flux 
implies cascade to the large scales. 

III. RESULTS 
A. A pilot run 

The first run in table U serves as a basic run to which 
all other runs are compared. For this reason this run is 
examined in more detail. Figure Q] shows the evolution of 
the kinetic and magnetic energy as a function of time. As 
can be seen the magnetic energy grows and saturates very 
fast at a relatively small amplitude. The kinetic energy 
on the other hand, after an initial fast growth transitions 
to a slower increasing phase. Up until the end of the 
numerical simulation this slow growth persists. The rea- 
son for this growth is the inverse cascade of the kinetic 
energy that accumulates energy in the large scales. 

This inverse cascade is demonstrated more clearly in 
figured This figure shows the parallel and perpendicular 
energy flux normalized by the total energy injection rate. 
The perpendicular energy flux is positive for wavenum- 
bcrs larger than the forcing wavenumbcr (direct cascade) 
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FIG. 3. Kinetic (top panel) and magnetic (bottom panel) 
energy spectra of run Rl. The solid lines correspond to the 
averaged spectra E u (k±), Et(k±) while the dashed lines indi- 
cate the zeroth component (fen = 0) of the two dimensional 
energy spectra E u (0, k±) and Eb(0,k±). The straight lines 
show for reference the power law spectra k~ 5 ^ 3 , k~ 2 , k~ 3 . 
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while for smaller wave numbers a negative constant en- 
ergy flux (inverse cascade) can be seen. On the other 
hand the parallel energy flux, shown by the dashed line, 
is small and always positive (ie direct). 

The presence of an inverse cascade can also indicated 
by looking at the energy spectra at late times. The top 
panel of figure [3] shows the kinetic energy spectrum E u 
of run Rl averaged over several outputs close to the end 
of the simulation. It can be clearly seen that most of the 
energy is concentrated in the large scales. The dashed 
line in this panel shows E u (0,k±). At large scales this 
line is identical to the E u spectrum thus the energy in 
these scales is mostly contained in the 2D-modes fcii = 0. 
This means that the flow in the large scales is almost 2D. 
(Here the flow is referred to as 2D in the sense that u has 
no dependence on the ^-direction and not that the u z 
component is absent.) On the other hand at the small 
scales E u (0, k±) is significantly smaller than E u thus the 
2D-modes contain only a small fraction of the energy and 
therefor the flow is three-dimensional. 

The bottom panel of fig [3] compares the magnetic en- 
ergy spectra Eb(0, k±) and E\>. Unlike the velocity field 
the magnetic field remains strongly three-dimensional for 
all scales since Eb(0,k±) <C Ej,. The amplitude of the 
magnetic energy is much smaller than that of the kinetic 
energy in the large scales but of the same order in the 
small scales. This is essential for the presence of the 2D- 
inverse cascade. If the magnetic field fluctuations were 
strong enough in the large scales the flow would behave 
as a 2D-MHD flow with a direct cascade. 

The fc -5 ' 3 scaling prediction for the 2D inverse cas- 
cade, the fc -3 for the direct 2D cascade and the fc -2 pre- 
diction of WTT are shown as a reference. The observed 
spectra are compatible with fc -5 / 3 in the large scales and 
k~ 2 in the small scales however the inertial ranges in the 
examined flow are too small to be conclusive. 



B. Guiding magnetic field strength 

As a next step the dependence of the inverse cascade, 
observed in Rl, on the amplitude of the uniform magnetic 
field is examined. Runs R2 and R3 have all parameters 
similar to run Rl but different value of the magnetic field 
amplitude. The flux of energy in both directions for the 
runs R1,R2 and R3 arc compared in figurc|4] As expected 
the amplitude of the uniform magnetic field has a drastic 
effect on the energy flux. The top panel of this figure 
shows II_i_(fcx)- R2 (dashed line) that has smaller value 
of V A than run Rl (solid line) has no inverse cascade and 
a stronger direct cascade. R3 (dashed-dot line) that has 
larger value of V A has on the contrary a stronger inverse 
cascade and a weaker forward cascade. The bottom panel 
of figured] shows the energy flux in the parallel direction. 
As the magnetic field is increased the flux to large k z is 
decreased. This expected since in the V A = oo limit there 
is cascade only in perpendicular direction. 

The spectra for these runs are compared in figure [S] 
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FIG. 4. Top panel: The energy flux in the perpendicular 
direction for Rl (V A = 5, solid line), R2 (V A = 2, dashed 
line), R3 (V A = 10, dashed-dot line). Bottom panel: The 
energy flux in the perpendicular direction for the same runs. 




10' 








10° 


R3 


,'"\ 




a 

UJ o 

10" 2 


_ R1 






LU 

io- 4 


- R2 




*\ VI 


10~ 6 




UK) 








10 


100 



FIG. 5. Top panel: The kinetic energy spectra of E u (k±) 
(solid line), and E u (0, k±) (dashed line) for R3, V A = 10 
(top lines), Rl, V A — 5 (middle lines), R2, V A = 2 (bottom 
lines). Bottom panel: The kinetic energy spectra E u (k±) 
(solid line), compared to the magnetic energy spectra Eb(k±) 
(dashed line) of the same runs and with the same order. The 
spectra have been shifted for reasons of clarity. 
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FIG. 6. Top panel: The energy flux in the perpendicular 
direction for Rl (8 < kf < 10, solid line), R4 (4 < k k < 5, 
dashed line), R5 (16 < kf < 20, dashed-dot line). Bottom 
panel: The energy flux in the perpendicular direction for the 
same runs. 



The top panel of this figure shows the kinetic energy 
spectra E u (k±) and E u (0,k±). The spectra have been 
shifted for reasons of clarity. In the two runs Rl and R3 
that showed an inverse cascade, energy is concentrated 
in the largest scales. What can also be observed is that 
as V A is increased the flow comes closer to a two dimen- 
sional flow. For R2 for which V A = 2 and no inverse 
cascade is observed the flow is far from two dimensional 
even at the largest scales. 

The bottom panel panel of figure [5] compares the ki- 
netic energy spectra E u (k±) (solid line) with the mag- 
netic energy spectra Eb(k±) (dashed line). As the uni- 
form magnetic field is increased the magnetic fluctuations 
are decreased compared to the velocity fluctuations. Note 
that for R3 the magnetic fluctuations are almost negli- 
gible in all scales while for R2 the two fluctuating fields 
are in equipartition. A possible interpretation of for this 
behavior is the following. Since in these runs there is 
no forcing for the magnetic field the magnetic fluctua- 
tions can be generated only by the stretching of field lines 
of the uniform component or by a dynamo mechanism. 
However since the flow comes close to a 2D flow as V A 
is increased neither of these mechanisms is possible. The 
dynamo mechanism however could depend on the mag- 
netic Reynolds number that for these runs it is relatively 
small. 



C. Forcing scale 

The second parameter we examine is the forcing scale 
kfL. This parameter is important because it controls the 
number of modes that satisfy the quasi-resonance condi- 
tions (HI),©- The energy flux of runs R4 with 4 < k f < 5 
and R5 with 16 < kf < 20 are compared to the energy 
flux of Rl with 8 < kf < 10 in figure [5] The top panel 
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FIG. 7. Top panel: The kinetic energy spectra of E u (k±) 
(solid line), and E u (0, k±) (dashed line) for R5, (top lines), 
Rl, (middle lines), R4, (bottom lines) Bottom panel: The 
kinetic energy spectra E u (k±) (solid line), compared to the 
magnetic energy spectra Eb(k±) (dashed line) of the same 
runs and with the same order. The spectra have been shifted 
for reasons of clarity. The arrows indicate the location of the 
forcing. 



again shows the energy flux in the perpendicular direc- 
tion while the bottom panel shows the energy flux in the 
parallel direction. 

All flows show an inverse cascade in the perpendicular 
direction. R4 has an inverse cascade of the same ampli- 
tude with Rl while R5 that is forced in smaller scales 
has a weaker inverse cascade. This is somehow expected 
since when the forcing is in smaller scales the system is 
closer in violating condition (J2J for 2D behavior. Note 
also that R4 has a larger flux in the parallel direction. 

The spectra for these runs are shown in figure [7J All 
runs have most of the kinetic energy concentrated in 
the large scales that behave like 2D-hydrodynamic flows: 
E u (kx) ^ E u (0,k±) (top panel) and E u (k±) > E b (k±) 
(bottom panel). The scales smaller than the forcing 
scale on the other hand behave like 3D-MHD flows with 
E u (k±) > E u (0, fc ± ) and E u (k±) ~ E b {k ± ). 
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FIG. 8. Top panel: The energy flux in the perpendicular 
direction for Rl (M = 0, solid line), R6 (M = 0.4, dashed 
line), R7 (M = 0.6, dashed-dot line). Bottom panel: The 
energy flux in the perpendicular direction for the same runs. 



D. Mechanichal and electro-motive forcing 

The inverse cascade observed in some of the discussed 
runs is a property of 2D hydrodynamic turbulence that 
is not present in 2D-MHD turbulence. The reason it 
appears in the previous runs is that the amplitude of 
the magnetic field fluctuations in the large scales re- 
mains weak. This effect could possibly be destroyed by 
a large scale dynamo at larger magnetic Reynolds num- 
bers. Leaving this possibility open the existence of the 
inverse cascade is investigated when magnetic field fluc- 
tuations are amplified by an electro-motive force. This is 
examined in runs R6 and R7 where both the mechanical 
and the electro-motive force are present. 

Figure E] shows the energy flux for runs R6 with M = 
0.4 and R7 with M = 0.6 compared with RI for which 
M = 0. The introduction of the electro-motive force 
(M ^ 0) destroys the inverse cascade in the perpendicu- 
lar direction (top panel) , while little change is observed in 
the parallel direction (bottom panel). This change indi- 
cates that the system transitions from an hydrodynamic 
2D state to a forward cascading MHD state. 

This is further confirmed by looking at the energy spec- 
tra in figure [HJ The top panel compares again the kinetic 
energy spectra E u and E u (0,k±). The excess of kinetic 
energy that is present in the large scales for run Rl, is 
absent in runs R6 and R7, verifying further the absence 
of the inverse cascade in the presence of an electro-motive 
force. Note that in all runs the large scales are still two- 
dimensional E u ~ -E^O, k±) (top panel) but the condition 
E u {k±) ^ &b(k±) is true only for Rl (bottom panel). 
This indicates that the absence of the inverse cascades 
for runs R6 and R7 is not because the flow stops behav- 
ing like a 2D flow, but rather because it starts behaving 
like a 2D-MHD flow. 

In 2D-MHD flows however there is an inverse cascade 




FIG. 9. Top panel: The kinetic energy spectra of E u (k±) 
(solid line), and E u (0,k±) (dashed line) for Rl, M = (top 
lines), R6, M = 0.4 (middle lines), R7 M = 0.6 (bottom lines) 
Bottom panel: The kinetic energy spectra E u (k±) (solid line), 
compared to the magnetic energy spectra Et(k±) (dashed 
line) of the same runs and with the same order. The spectra 
have been shifted for reasons of clarity. 
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FIG. 10. The vector potential spectrum E A for three different 
times. 
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FIG. 11. Top panel: The energy flux in the perpendicular 
direction for Rl (solid line), R8 (dashed line), R9 (dashed- 
dot line). Bottom panel: The energy flux in the perpendicular 
direction for the same runs. 



of the square of the vector potential that is a conserved 
quantity. If the flow in the runs R6 and R7 behave like a 
2D MHD flow in the large scales such a cascade should be 
observed. However a flux for the squared vector potential 
in three dimensions can not be uniquely defined since it is 
not a conserved quantity. Nonetheless, we plot the vector 
potential spectra for three different times from run R7 
in figure [TUJ The vector potential a is defined so that 
b = V x a and V ■ a =0. Its spectrum is then defined as 



(3) 



where the sum is restricted in the wavenumbers k± < 

+ ky < fcj_ + 1 and a^ is the Fourier transform of 

a. It can be seen that as time progresses the vector po- 
tential is increasing in the large scales. It is noted that 
a quasi-conservation of the square of the vector potential 
has been observed in [3Cj [Hi for three-dimensional ideal 
reduced-MHD. 



E. Isotropic and anisotropic forcing 

The last parameter that we vary is the isotropy of the 
forcing. Unlike the previously examined runs for which 
all Fourier modes within a spherical shell arc uniformly 
forced in runs R8 and R9 the modes fc|| = and k± = 
respectably are forced preferentially. In particular the 
amplitude of a Fourier mode Fk of the forcing inside the 
chosen spherical shell was proportional to: 



R8 : |F k | oc 
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R9: |F k | oc ||, (4) 
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where k 2 = k 2 x + k 2 + k 2 z . Thus in run R8 the 2D modes 
(k z = 0) are forced preferentially while in run R9 these 



FIG. 12. Top panel: The kinetic energy spectra of E u (k±) 
(solid line), and E u (0,k±) (dashed line) for R8, (top lines), 
Rl, (middle lines), R9 (bottom lines) Bottom panel: The 
kinetic energy spectra E u (k±) (solid line), compared to the 
magnetic energy spectra Eb(k±) (dashed line) of the same 
runs and with the same order. The spectra have been shifted 
for reasons of clarity. 



2D modes are not forced at all and most forcing is at the 
high k z modes. 

The energy flux for these runs is shown in figure 1111 
The energy flux of run Rl and R8 are very close for both 
directions. The reason for the slightly smaller flux of 
R8 at the largest scales is because R8 was evolved for a 
shorter time than Rl. R9 however has no inverse cas- 
cade in the perpendicular direction (top panel). This is 
not surprising since for this run the 2D modes are not 
forced at all. In the parallel direction there is also a 
drastic change. R9 shows an inverse cascade from the 
large k z wavenumbers to the 2D k z — modes. It is 
noted that the flux towards the large parallel scales was 
strongly fluctuating taking both positive and negative 
values. Only after averaging several files the result shown 
in figure [TT] was obtained. 

The spectra for these runs are shown in figure 1121 
Again not a lot of difference can be seen between run 
Rl and R8. Both are close to two-dimensional in the 
large scales and three dimensional in the small scales (top 
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panel), and both have weak magnetic energy in the large 
scales but are close to cquipartition in the small scales 
(bottom panel). R9 on the other hand is three dimen- 
sional for all scales and kinetic and magnetic energy are 
almost identically equal at all scales. This last remark in- 
dicates that Alfven-waves (for which u = ±b) dominate 
the turbulence. 



IV. SUMMARY AND DISCUSSION 

In this work we have shown that under certain condi- 
tions an MHD-flow in the presence of a strong magnetic 
field can behave like a two dimensional flow in the large 
scales while like three dimensional (possibly weak tur- 
bulence) in the small scales, much like strongly rotating 
fluids. In the large scales it was found that the magnetic 
fluctuations were suppressed and a 2D inverse cascade of 
energy developed with energy accumulating in the large 
scales. 

This inverse cascade however is sensitive to various pa- 
rameters. If the uniform magnetic field amplitude is de- 
creased sufficiently, or if the domain size is increased (or 
equivalently the forcing scale is decreased) so that the 
condition © for two-dimcnsionalization no longer holds 
the flow recovers its 3D behavior and cascades the energy 
forward. Furthermore, in the case that an electro-motive 
force is introduced although the flow remains 2D in the 
large-scales the magnetic fluctuations are no longer sup- 
pressed and the flow behaves like a 2D-MHD fluid with 
a direct energy cascade. At the same time an indication 
of an inverse cascade of the squared vector potential was 
observed. Finally, absence of an inverse cascade in the 
perpendicular direction was also observed when the sys- 
tem was forced only in the large fcy modes. In this case an 
inverse cascade in the direction parallel to the magnetic 
field was observed. 

All the simulations presented here that exhibited an 
inverse cascade, were stopped before the largest scale of 



the system was reached due the computational cost. If 
the runs were continued for a longer time as the energy 
and the perpendicular wave number are increased it is 
possible that a point will be reached that the strength 
of the magnetic field will not be sufficient to stop three 
dimensional instabilities from breaking the two dimen- 
sional constrain. Such a transition point is expected to 
appear when the eddy turn over frequency Ukk± becomes 
smaller than the Alfven frequency B^k^ . In the inertia 

range of the inverse 2D cascade where the k ± 5 ^ 3 scaling 
is expected the eddy turn over frequency is decreasing as 
larger scales are reached. On the contrary the smallest 
Alfven frequency B^k^ remains independent of k±. Thus 
the ratio Ukk_\_/ B$k\\ will decrease as the cascade pro- 
ceeds and it is not expected that such a transition point 
will exist in the inertial range, instead as the cascade 
proceeds the flow will come closer to a 2D flow. How- 
ever since there is no large-scale damping mechanism to 
dissipate the energy when the largest scale of the sys- 
tem is reached the energy would pile up in this scale. In 
this case the eddy turn over frequency will increase and 
eventually the criterion @ for two-dimensionality will be 
violated. Then energy could possibly return to the small 
scales as weak or strong turbulence. Similar scenario for 
the fate of the inverse cascade of rotating turbulence has 
been proposed in [32| . Such scenarios however need to 
be verified by numerical simulations and experiments. 
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